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In this work, we theoretically examine recent 
pump/probe photoemission experiments on the strongly 
correlated charge-density-wave insulator TaS2- We 
describe the general nonequilibrium many-body for- 
mulation of time-resolved photoemission in the sudden 
approximation, and then solve the problem using 
dynamical mean-field theory with the numerical renor- 
malization group and a bare density of states calculated 
from density functional theory including the charge- 
density-wave distortion of the ion cores and spin-orbit 
coupling. 



We find a number of interesting results: (i) the bare band 
structure actually has more dispersion in the perpendicu- 
lar direction than in the two-dimensional planes; (ii) the 
DMFT approach can produce upper and lower Hubbard 
bands that resemble those in the experiment, but the up- 
per bands will overlap in energy with other higher en- 
ergy bands; (iii) the effect of the finite width of the probe 
pulse is minimal on the shape of the photoemission spec- 
tra; and (iv) the quasiequilibrium approximation does not 
fully describe the behavior in this system. 



1 Introduction Strongly correlated electronic systems 
are some of the most interesting and exotic systems in the 
quantum world. They are compounds where the outermost 
electrons interact so strongly with each other that nonin- 
teracting band theory no longer applies, and they can ex- 
hibit novel quantum-mechanical effects. Most experimen- 
tal and theoretical techniques for strongly correlated elec- 
trons have been applied in equilibrium, or in the linear- 
response regime, where we have a wealth of knowledge 
and understanding about how these systems behave. Strongly 
correlated electron systems that are driven into a nonequi- 
librium state are much less understood, and represent a new 
frontier for the condensed-matter physics community. One 
of the most important experimental techniques that allows 
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access to the nonequilibrium regime are time-resolved pump/ 
probe studies, where an intense and ultrashort pulse of light 
hits the material, causing it to become excited into a nonequi- 
librium state, and then a much less intense but still ultra- 
short pulse is used to probe the system at different time 
delays after the pump pulse, to examine how the nonequi- 
librium dynamics lead the system to relax back to equilib- 
rium. In this work, we will focus on time-resolved pump/ 
probe photoemission experiments that have recently been 
completed in TaS 2 |ffl. 

The material TaS2 is an interesting strongly correlated 
commensurate charge-density-wave (CDW) insulator be- 
low an incommensurate/commensurate transition tempera- 
ture of about 180 K [2|. Conventional band-structure cal- 
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Figure 1 (Color online.) Structure of the CDW-ordered Ta 
plane in the IT3 phase at low temperature. There are three 
inequivalent Ta sites in the ordered phase, as indicated by 
the different colors. 



culations predict that the ordered phase is metallic, with a 
single band crossing the Fermi level, but experiment shows 
that the material is actually an insulator at low T, and it 
is believed the insulating state arises from Mott-Hubbard 
physics in this single band O (which was originally pro- 
posed to be close to the critical U for the Mott transi- 
tion). The CDW structure is actually quite complex in this 
system. Starting at temperatures above 1150 K, we find a 
(disordered) IT structure is the stable structure for TaS2; 
this structure is made up of S-Ta-S trilayer units, in which 
atoms in each layer sit on a triangular lattice [4j. Since the 
separation between S-Ta-S sandwiches is larger than the 
interlayer spacing within a sandwich, the structure is quasi 
two dimensional. As the temperature is lowered, TaS2 un- 
dergoes a series of phase transitions to structures that can 
be viewed as distortions of the basic IT structure. Of inter- 
est here is the low-temperature commensurate CDW phase 
(IT3), which appears below about 180 K. In the plane, this 
phase is characterized by a x \/l3 supercell in which 
the Ta atoms group into clusters of 13 sites in a 'Star-of- 
David' arrangement, as shown in Fig. Q] Some groups re- 
port that the stacking sequence of S-Ta-S sandwich units 
has a period of 13 (4), while others suggest that the stack- 
ing is disordered J5). For this work, we employ a \/l3 x 
\/l3 x 1 supercell to model the IT3 phase; i. e., we assume 
uniform stacking in the z-direction. 

2 Many-body formalism for TR-PES The theoret- 
ical description for time-resolved (TR) pump/probe pho- 
toemission spectroscopy (PES) is complicated. For contin- 
uous beam photoemission, the signal is proportional to the 
so-called lesser Green's function (summed over k z if the 



system is not quasi two-dimensional) in frequency space, 
when one neglects the specific details of the matrix ele- 
ments involved in the coupling of the electrons to the light 
and in the sudden approximation (see below for details). 
For the nonequilibrium case, the formulation of the photoe- 
mission theory is more complicated. One first must evolve 
the system in the presence of the pump pulse to create 
the 'initial' nonequilibrium state on which the probe pulse 
acts. This is achieved by evolving the equilibrium ensem- 
ble of many-body eigenstates of the system Hamil- 
tonian H to the 'initial' ensemble of states {^(^o))} = 
{U(to, — oo)!^)}, where U is the unitary time-evolution 
operator in the presence of the pump pulse, determined 
by a time dependent Hamiltonian 7i pU mp(i) whose pre- 
cise form depends on the way one models the interaction 
of the pump radiation [represented by the the vector poten- 
tial A pump (r, t), with the t dependence including the turn- 
ing on and off of the pump,] with the electrons of the sys- 
tem 16) . The turning on of the probe pulse at time to sim- 
ilarly modifies H by adding in the Hamiltonian 7Y pro be(i) 
which includes the interaction of the vector potential of 
the probe radiation with the system. At some later time 
t, (typically after the probe pulse has been turned off) the 
initial nonequilibrium ensemble of states {|^(io)}} has 
evolved into the 'final' ensemble of states {\&„ (*))} = 
{J7(i,trj)|^n(£o)}}> with U(t,t ) being the unitary time- 
evolution operator that includes the presence of the probe 
pulse. The probability that a photoelectron with momen- 
tum k e = k e k e (in a momentum interval dk e and solid 
angle dflt ) is detected when the system is describable by 
this final nonequilibrium ensemble is given by 



lim 



(k e ) dk e dQj, 



(2tt) 3 



^P(t); P(t) = J2Pn\{^m;-keWn(t))\ 



(1) 

Here, as appropriate to photoemission from an experimen- 
tal sample with a surface, |S^ m ;k e ) is well approximated 
as a direct product of the many-body eigenstate \*P m ) of 
the initial (and final) time-independent system Hamilto- 
nian TL, and a high energy, one-electron scattering state 
of the system which propagates as a free electron state of 
momentum Kk e outside the sample. The eigenstate YP m ) 
in which the system is left is not determined in the experi- 
ment, and the initial state can be any one of the ensemble of 
initial states with probability p n = Z^ 1 cxp[— E n / (fcgT)] 
where E n are the corresponding energy eigenvalues, and 
Z = J2 n ex P[ — En/ (fcsT)] is the partition function. Hence 
Eq. (HJ includes an unconstrained sum over m, and a sum 
over n weighted by p n . This expression implicitly assumes 
the sudden approximation, where the photo-excited elec- 
tron rapidly moves out of the sample. 

We employ a probe shape function, s(t'), to capture 
the time dependence of the temporal profile of the probe 
pulse, including its turning on and off. M q (fc z , fe ez ;k||), 
the matrix element for the absorption of a photon of wave 
vector q (frequency w q = cq) and the ejection of an elec- 
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tron of wave vector k = (kii,fc z ) inside the system as a 
photo-electron of wave vector k e = (k||, k ez ) depends on 
the details of the modeling of the sample, especially its sur- 
face . To leading order in the perturbation Ti pro b e it ) due 
to the (weak) probe pulse (and using the factorization of 
l^mjkg)), we find 



|(^ m ;k e K(i))| 



J dfc z M q (fc z ,A: ez ;k|| 



x / dt' s{t')e- iut ' 



(2) 



x (* m \U\t',t )c k U(t',t )\*l(t )) 



with U being the exact unitary time evolution operator in 
the presence of the pump pulse, which therefore includes 
all its nonequilibrium effects. Here hw = 7kJ q — 
(hk e ) 2 1 (2m e ) — W is the energy of the excitation left in 
the system after the photoemission process and W is the 
work function. For probe photons of fixed direction and en- 
ergy, and for a given material, specifying ui determines k e , 
hence the probability P(t) in Eq. (JTJ is a function only of 
t, ui and k e . Using the properties of the time development 
operator, and the completeness of {I^Vn)} it is straightfor- 
ward to show that 



P{t,UJ, fc e ) 



with 



1 



dk' z / dfc z M q (fc z ,fc ez ;k||) 



x A^(fc z ,/c ez ;k||)/(t,w,fc e ;fe,fc z ) (3) 



I{t,ui,k e ;k z ,k' z ) = -i dt" tft / *(t")«(t / )e M *"~*' ) 
xG<(t',t"). (4) 



it 1 — t"). Furthermore, in a highly anisotropic layered sys- 
tem, its k z dependence can be neglected. Then the angle- 
resolved PES transition rate (transition probability per unit 
time, which is what is relevant as the continuous beam 
probe pulse photoemits electrons at all times at a constant 
rate) is proportional to 



lim lim 

t — >OC tQ — > — c 



I(t,u, k e ] 
(f - to) 



Here k' = (k||, k' z ) , and G < is the well known two-time 
(nonequilibrium) lesser Green's function (8| given by 

G< w {t',t") = iY;p n {# n \u(-ooX)<i>WJ) 

n 

x Ck^t'.-oo)!^) (5) 

= iZ- 1 Hr[e- w /( fc * r )4,(t")c k (t')] (6) 

where c^,(t") and Ck(t') are the electron creation and de- 
struction operators in the Heisenberg picture appropriate to 

^pump(^)' 

4,(t") = [/(-oo,i")4,[/(i",-oo); (7) 
Ck(i') = U{-oo, t')c k U(t', -oo). (8) 

In the nonequilibrium (pumped) case, this needs to be cal- 
culated using nonequilibrium Keldysh contour-ordered Green 
function techniques [9 | on the Kadanoff-Baym- Keldysh con- 
tour in the complex time plane. 

In case of continuous probe beam PES on a system in 
equilibrium, s(t) = 1 and G^(t', t") is only a function of 



(9) 

which is the standard result [with (w) being the spec- 
tral function and f{u>) — l/{l + cxp(ui/kBT)} the Fermi- 
Dirac distribution function] . 

We will also invoke the quasiequilibrium approxima- 
tion in this work. This approximation assumes that the elec- 
tronic system rapidly thermalizes after the pump pulse is 
turned off, but, because the time scale for thermalization 
with the phonons is much longer, the electronic system re- 
mains hot, and cools slowly during the course of the TR ex- 
periment JTJ . In this case, the PES is described by the equi- 
librium lesser Green's function, but at an effective elec- 
tronic temperature that depends on the time delay for the 
turning on of the probe pulse, convoluted with the probe 
pulse shape functions s(t) as in Eq. [4] Finally, we take 
the matrix elements M to be constants, and evaluate the 
Green's functions for a bulk system, rather than a system 
with a surface, in which case we have full translational in- 
variance, and k z = k z . This last approximation is perhaps 
the most serious one that we make, but it also is commonly 
done in theoretical treatments of PES. 

3 Numerical techniques Because we need to describe 
Mott-Hubbard physics we will use a dynamical mean-field 
theory (DMFT) approach. Since there is only one band at 
the Fermi level, this problem can be formulated as a sim- 
ple single band model with a complicated noninteracting 
DOS. (This tacitly assumes that any other higher bands re- 
main higher in energy than the upper Hubbard band that 
comes from the original noninteracting band at the Fermi 
energy, a result that may well not be true in TaS2 - see com- 
ments below.) That DOS is determined via band structure 
calculations that include both the CDW distortions of the 
ion cores and the spin-orbit coupling. 

The band structure calculations were done using the 
density-functional code VASP 11101 . a plane-wave-based all- 
electron code in which the electron-core interaction is treated 
using the projector augmented wave method ifTTl . The gen- 
eralized gradient approximation (GGA) was used for elec- 
tron exchange and correlation lfT2l . For self-consistent-field 
calculations and structural relaxations, the Brillouin zone 
was sampled using a 6 x 6 x 12 grid of k-points. Unless ex- 
s plicitly noted, calculations were carried out using the scalar 
relativistic Hamiltonian. 

The GGA optimized in-plane lattice constant of a = 
12.2 A is about 0.5% larger than the experimental value. 
For the out-of -plane lattice constant, the experimental c/a 
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Figure 2 (Color online.) Band structure of CDW-ordered 
TaS2, including the effects of the spin-orbit interaction. 
The zero of energy is set to the Fermi level. To facilitate 
comparison with results in Ref. iTPTl which ignored spin- 
orbit coupling, the bands are plotted along high-symmetry 
directions in the Brillouin zone of the undistorted IT struc- 
ture. The half-filled band (red) that crosses the Fermi level 
is strongly localized in the in-plane directions, and is dis- 
persive along k z (r to A). 



was used to mitigate effects arising from the simplified 
stacking sequence assumed. Atomic positions were relaxed 
in the scalar relativistic approximation, holding the lattice 
constants fixed. The b- and c-site Ta atoms relax towards 
the center a site, with the a-b and a-c distances contract- 
ing by 5.3% and 3.8% from their respective values in the 
undistorted triangular lattice. These results are close to the 
experimentally measured distortions. For the S atoms, the 
relaxations were found to be larger in the z direction than 
within the planes. The S planes develop a slight pucker, 
with S atoms close to the central Ta a sites displaced out- 
wards in the z direction. 

In the undistorted IT structure, a previous density func- 
tional calculation has shown that the electronic bands have 
a two-dimensional character in which the Ta d bands cross- 
ing the Fermi level are only weakly dispersing in the k z 
direction, consistent with the quasi 2D nature of the lattice 
[ 1 3 1 . That study also found that the formation of the Star- 
of-David clusters causes significant changes in the bands 
near the Fermi level. Indeed, in our calculation, the CDW 
distortion causes the Ta d band that crosses the Fermi level 
to become essentially one-dimensional, with a width of 
about 0.5 eV arising almost exclusively from dispersion 
in the k z direction (see Fig. 2). This is a nonbonding d z i 
band with primary weight on the central Ta a site. The Ta 
clustering causes this state to become localized in the in- 
plane directions, which suggests that the electronic struc- 
ture should transition from being two-dimensional to zero- 
dimensional. However, the outward bulging of the S atoms 
around the central Ta a sites that accompanies the forma- 
tion of clusters enhances the interactions between tri-layer 
units, resulting in a one-dimensional band instead. 



Recently, it has been suggested that in addition to the 
CDW distortion, the spin-orbit effect could have a large 
influence on the bands near the Fermi level in TaS2- A 
2D tight-binding model found that the spin-orbit interac- 
tion causes the nonbonding Ta d band to split off from 
a manifold of higher-lying bands |fl4l . It becomes a very 
narrow (less than 0.1 eV in width) half-filled band right 
at the Fermi level, primarily of d xy , d x 2_ y 2 character on 
the Ta a site. We have carried out density-functional cal- 
culations for the 3D system, treating the spin-orbit effect 
self-consistently. The results are shown in Fig.|2] As in the 
tight-binding results, we find that the spin-orbit interaction 
splits off the a-site nonbonding band from higher-energy 
bands. However it still has a width of about 0.4 eV due to 
k z dispersion. Its character remains mainly d z 2, although 
the d xy and d x 2_ y 2 weight is significant at the top of the 
band. 

The density of states (DOS) of the split-off one-dimensional 
band is shown in Fig. [3] with the {7 = label. Its double- 
peaked shape can be understood from the k z dispersion: 
near the bottom of the band, the dispersion is quadratic, 
leading to a DOS that goes roughly like (E— Eq)~ x / 2 \ near 
the middle of the band, the dispersion is roughly linear, giv- 
ing a nearly constant DOS; and at the top of the band, the 
band is relatively flat in all three directions, yielding a peak 
in the DOS. 

In spite of the fact that the dispersion of the CDW- 
distorted bands (with spin-orbit coupling) have a quasi one- 
dimensional character to them, we nevertheless apply the 
DMFT approach to determine the Mott-Hubbard physics, 
even if the results will only be approximate. The DMFT 
calculations are now straightforward given the noninter- 
acting DOS. We make the simplifying assumption that the 
Coulomb interaction will be between the two different (de- 
generate) spin-orbit states at a given local site and use the 
DFT band structure DOS (for the single band that crosses 
the Fermi level) in the Hilbert transform that determines 
the local Green's function on the lattice from the momentum- 
dependent Green's function with a local self-energy. Then, 
the conventional numerical renormalization group (NRG) 
impurity solver can be used to solve for the Green's func- 
tion of the effective impurity problem with the given value 
of the interaction U. The NRG approach we employ is 
the standard algorithm ifTBI — we take A = 1.8 and keep 
800 states per iteration on the Wilson chain. The DMFT 
equations for the retarded Green's function are solved via 
the self-consistent iterative approach of Jarrell |16|. The 
PES signal is found by multiplying the spectral function 
by —2if{w) to get the lesser Green's function, then Fourier 
transforming to real time to find G < (t — t'). Next we eval- 
uate Eq. (@J, which is proportional to the TR-PES spectra 
given the above assumptions. We take the shape function 
s(t) to be a Gaussian s(t) = exp[-(i - i) 2 / T 2 } / YT \fii) 
with a width J7 equal to 80 fs and a probe delay time set at 
t (we assume the center of the pump pulse occurs at t = 0). 
In the quasiequilibrium approximation, the lesser Green's 
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Figure 3 (Color online.) Local many-body density of 
states for the single-band Hubbard model (at T = 50 K) 
with the noninteracting DOS given by the band structure 
calculations described in the text. The metal-insulator tran- 
sition takes place around U f=a 0.55 eV, and the best param- 
eter value for TaS2 is U ~ 0.7 eV. The dashed black lines 
show the upper and lower bands in the band structure that 
lie above and below the band that crosses the Fermi en- 
ergy; note how the upper and lower Hubbard bands merge 
with them more and more as U increases (even the nonin- 
teracting band is overlapping with neighboring bands near 
the lower band edge). 



function depends only on the time difference, so the convo- 
lution theorem can be used, and the final PES signal can be 
determined directly by a single integral in frequency space 



dvG< 



(^-^))| S »| 2 /27r, 



(10) 



where s(y) is the Fourier transform of s(t'). Since a pulse 
of 80 fs has a width on the order of 8 meV, this convolu- 
tion will not modify the smooth structure in the DOS or 
the PES which vary over a frequency scale much larger 
than 0.01 eV. Hence we neglect the "windowing effect" of 
the probe envelope function s(t) in our analysis below. We 
fix our calculations to half-filling and vary the interaction 
strength to tune it for TaS2, and then only vary the tempera- 
ture. Our PES signals are normalized so that the integrated 
weight in all spectra are identical. 

We begin by calculating the local DOS at T = 50 K 
to set the value of U. The results are shown in Fig. [5] 
Note how the metal-insulator transition takes place near 
U = 0.55 eV. We "fit" to the experimental data by pick- 
ing the distance between the peaks of the upper and lower 



Hubbard bands to be approximately 0.6 eV, as seen in the 
TR-PES studies 1 1 . j. This result also produces a gap on the 
order of 0.125-0.15 eV, consistent with optical conductiv- 
ity experiments |[T7l . In addition, one should note that the 
upper and lower Hubbard bands are merging with the next 
higher and next lower bands in the band structure, indicat- 
ing that the single-band model may be inadequate for this 
system (particularly for the upper Hubbard band). 

4 Results The experimental continuous beam (equi- 
librium) PES studies on TaS2 in the low-T CDW phase 
show interesting behavior fl8l . The angle-resolved PES 
shows little dispersion of the lower Hubbard band peak as a 
function of momentum, consistent with the dispersion ly- 
ing predominantly in the k z direction (the k z dispersion 
is seen in PES studies that vary the photon energy lfl3ll ). 
Furthermore, as one approaches the chemical potential, the 
angle-resolved PES curves tail off toward zero with an al- 
most linear dependence on frequency. This does not dis- 
play the expected band gap in the insulating phase, where 
the spectra should go to zero before one hits the chemical 
potential. It is argued in the experiment that this is occur- 
ring because the system has defects which pin the chemi- 
cal potential to the upper edge of the lower Hubbard band. 
If so, then there would be no way to see the gap in any 
continuous-beam (equilibrium) PES experiment, because 
one cannot reach the high temperatures needed before the 
CDW order disappears. In TR-PES experiments, one can 
access the higher bands by pumping energy into the sys- 
tem with the pump pulse which can excite electrons across 
the CDW gap and then be imaged. But the TR-PES spectra 
also do not show a clear signature of a gap in the data, as 
we illustrate below. 

The TR-PES measurements on TaS2 work with low- 
energy photons HI, and hence the electronic excitations 
are limited to lie in a small volume about the zone-center 
in the Brillouin zone. It is argued, if the system is quasi 
two-dimensional, that one would then examine the k = 
spectral function. We will show some results with just such 
an approach, although, the quasi one-dimensional disper- 
sion along k z , indicates that one should average instead 
over k z , which is approximately the local Green's function 
when the dispersion is quasi-one-dimensional. So we also 
examine the local Green's function with our simplified ap- 
proach. 

Since we have set the Hubbard model interaction at 
U = 0.7, we can determine the TR-PES spectra by simply 
calculating the local (or k = 0) equilibrium lesser func- 
tion at different temperatures and plotting as a function of 
frequency. Different effective temperatures correspond to 
different time delays, with the hottest temperatures being 
the shortest time delays. The results of just such a calcula- 
tion are shown for the local Green's function in Fig. @] and 
for the k = Green's function in Fig. [5] Here we plot the 
lesser function for different quasiequilibrium temperatures. 
The chemical potential is fixed by the requirement that the 
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Figure 4 (Color online.) Time-resolved photoemission 
spectroscopy in the single-band quasiequilibrium approx- 
imation for TaS2- Here we plot the local lesser Green's 
function to represent the sum over all k z values with the 
quasi one-dimensional dispersion. Note how one always 
sees a strong insulator, and that one can image the upper 
Hubbard band as the temperature becomes on the order of 
half the bandwith and above. Note further that the 50 K 
curve and the 175 K curve are indistinguishable. 



Figure 5 (Color online.) Time-resolved photoemission 
spectroscopy in the single-band quasiequilibrium approx- 
imation for TaS2- Here we plot the k = lesser Green's 
function. In this case, the relative weight of the upper and 
lower band features are closer to experiment than the local 
results and the tailing off of the signal as we approach the 
chemical potential from negative energies is also similar to 
what is seen in experiment, but the peak location is pushed 
much lower in energy and sits near —0.5 eV rather than 
—0.2 eV, as seen in experiment. 



filling remain exactly one in the band. The chemical po- 
tential shifts in such a way to push the peaks to the left 
in the figures for higher temperature. Such a shift is not 
seen in the experimental data. If we instead assume that 
the chemical potential is pinned by defect states, then we 
would shift the curves back to the right (by hand) so they 
all meet at a characteristic frequency. This is the behavior 
seen in experiment and could be reproduced in the theory if 
we allow ourselves to shift the frequency axis for each case 
or if we include defect states in our analysis. In addition to 
the movement of the peaks with the changing chemical po- 
tential, we also clearly see insulating behavior for small T, 
and we always see a double peak structure when we can 
image the lower and upper Hubbard bands separately at 
high T. Note that the upper Hubbard band is always small 
in weight compared to the lower band due to the Fermi 
factor. The experimental TR-PES data does not show such 
behavior. It instead has an almost monotonically decreas- 
ing signal as one moves across the chemical potential and 
beyond, with no real gap structure to be seen anywhere. 
This behavior is difficult to reconcile with a single-band 
Mott-Hubbard type of picture. The behavior for some of 
these features is reproduced better in the local approxima- 
tion than the k = approximation (like the location of the 



peak height), while other behavior is produced better in the 
k = approximation (like the tailing off of the signal as 
one approaches the chemical potential, the lower relative 
weight of the upper Hubbard band, and the difficulty in di- 
rectly imaging the peak in the upper Hubbard band). 

We have a couple of conjectures about what is happen- 
ing in this system, but we do not yet have definitive an- 
swers to these questions. First off, the single-band model 
approximation is probably the most serious one that we 
made in this work. This is because the higher-energy bands 
in the band structure lie just slightly above the noninter- 
acting DOS upper band edge. So the upper Hubbard band, 
which is pushed upwards in energy relative to the noninter- 
acting band edge, is expected to lie in or above those first 
lower bands (the situation is not quite as bad for the k = 
spectra, because the k = band energies are pushed a bit 
higher). This implies the system looks more like a charge- 
transfer insulator than a Mott-Hubbard insulator, and the 
upper bands in the band structure must be taken into ac- 
count in the analysis. If we instead ignore that problem, 
then we still need to explain why we do not see a clear sig- 
nature of the upper Hubbard band in the data. Here there 
are three other effects that can be playing a role. One could 
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have a situation where the nonequilibrium behavior is not 
described well by the quasiequilibrium approximation, and 
hence the "filling in" of the gap is a nonequilibrium ef- 
fect. Another possible explanation could arise from surface 
states which may lie at slightly different energies and could 
enhance the signal in the gap, making it look like there is 
no gap, when there actually is one in the bulk. Finally, the 
energy and momentum dependence of the matrix elements 
could be playing a role. When one compares the TR-PES 
studies with the higher-energy continuous beam studies, 
one can see a difference in the relative peak heights of the 
lower Hubbard band and of the even lower energy filled 
bands, which must arise from matrix-element effects. 

5 Conclusions In this work, we have examined a the- 
oretical description of the time-resolved photoemission spec- 
troscopy in the strongly correlated CDW insulator lT-TaS2 
based on a simplified single-band model that employs DFT 
+ DMRG(NRG) techniques. The model can be solved di- 
rectly, and does display some of the qualitative features 
of the experiments, but misses a number of key elements: 
First, the experiments never show a true insulator, while 
the theory shows well defined insulating behavior even at 
high temperature, where a lower and upper Hubbard band 
can still be seen. Second, this approximation is neglect- 
ing the fact that the upper Hubbard band lies at energies 
that are typically above the lower band edge of the higher 
energy unoccupied band states of the system. When this 
occurs one really needs to consider a multiband model. 
Third, we find the DFT-based bands, with the CDW dis- 
tortion and the spin-orbit coupling included, display quasi 
one-dimensional behavior perpendicular to the planes, and 
essentially localized behavior within the planes, indicating 
that averaging over k z is always necessary (even though 
many features in the experiment are represented better by 
the k = results with no k z averaging). It is also possi- 
ble that one need to take into account nonequilibrium ef- 
fects beyond the quasiequilibrium approximation (which 
may fill the gap in the DOS faster), include the momentum 
or energy dependence to the matrix elements rather than 
treating them as constants, and incorporate surface effects, 
in particular the effects of surface states which could also 
mask the insulating behavior. Clearly further work is re- 
quired to start to resolve and better understand lT-TaS2- 
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